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We have studied the properties of clusters of colloidal magnetic particles generated from a 2D 
aggregation model with dipolar interparticle interactions. Particles diffuse off-lattice, experiencing 
dipolar interations with the already attached particles until either they stick to the cluster or wander 
far away and are removed. Our results are interpreted in terms of a fractal dimension that is a 
monotonically decreasing function of the temperature, varying between a definite value close to 1 
at T = 0, and the limit T ^ oo, corresponding to free diffusion-limited aggregation. By analyzing 
orientational correlation functions, an ordered state is found at low temperatures; this state is 
destroyed by the fractal disorder generated at high T. Our study could be relevant in understanding 
aggregation of dipolar colloids and phase transitions in Langmuir monolayers. 
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I. INTRODUCTION 

Fractal growth phenomena [|l| have been a particularly 
active field of physics in the last decade, due to the poten- 
tial applications to many disciplines, in particular to the 
physics of colloids. In the context of computer models, 
the most noteworthy are the diffusion-limited aggrega- 
tion (DLA) model |^ and the cluster-cluster aggregation 
model [H , which indeed can describe the fractal structure 
of colloidal aggregates 

Most often the above mentioned models assume very 
short range interparticle interactions; usually a hard-core 
potential plus an infinite well on the surface of the par- 
ticles. They are therefore appropriate to describe aggre- 
gation with interactions strongly decaying with distance. 
However they fail to represent aggregation in the pres- 
ence of long-range forces. Some examples of such kind 
of processes are the aggregation of particles subject to 
dipolar forces. Actual experiments have been conducted 
in ferrofluids and the so-called magnetic holes ||] ; part of 
their interest lies in the extreme simplicity of the exper- 
imental setup, which can be easily performed with mag- 
netized microspheres |^]. On the other hand, there is a 
considerable interest in the study of phase transitions in 
Langmuir monolayers Q . Far from equilibrium, the con- 
densed phase grows at the expense of the liquid phase, 
forming clusters of different shapes. The phospholipids 
constituting the Langmuir monolayer experience repul- 
sive dipolar interactions, which must play a major role 
in determining the morphology of the condensed aggre- 
gates. 

Some authors have extended the classic models in order 
to take interparticle interactions into account. The mod- 
iflcations proposed so far are either deterministic 
or random pG-nj, using Monte Carlo or Langevin dy- 



namics methods. However, because the implementation 
of long range interactions is extremely time-consuming, 
only modest-sized aggregates can be grown (up to 128 
particles in |^,|lj,|l5[ ) . Therefore the results do not al- 
low categorical conclusions to be made about the fractal 
properties of the dipolar clusters. 

Our purpose in this paper is to extend the particle- 
cluster aggregation model in order to include fully 
anisotropic attractive dipolar interparticle interactions. 
Our algorithm is considerably fast, allowing us to gen- 
erate clusters up to 10000 particles at zero temperature, 
in a reasonably short time. In Section II we describe the 
technical details of the algorithm. Section HI discusses 
the evolution as a function of temperature of the fractal 
dimension and the order induced on the orientation of the 
dipoles by their reciprocal interactions. Our conclusions 
are discussed flnally in Section IV. 



II. CLUSTER FORMATION ALGORITHM 

We consider the two-dimensional aggregation process 
of magnetic particles of diameter d and magnetic moment 
jl = jiu^ with ji being the magnetic moment strength 
and u a unit vector oriented along its direction. The 
dipolar energy between two particles i and j, located at 
the positions fi and r^, respectively, is £ij = /J.'^Eij, Eij 



being the dimensionless dipolar energy 

Eij = {ui-Uj- 3{ui ■ fij)(uj ■ fij)/r^j} /rl 



(1) 



and fij = fi — fj. 

The simulation starts with a seed particle located at 
the origin of coordinates, bearing a randomly oriented 
tridimensional vector ui^ parallel to the plane of growth. 
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The following particles are released from a random posi- 
tion on a circle of radius Rin centered on the seed. The 
particles have assigned a vector Ui oriented at random. 
Each particle undergoes a random walk until it either 
contacts the cluster or moves away from the origin a dis- 
tance greater than i^out- In this case, the particle is re- 
moved and a new one is released from the circle surround- 
ing the seed. We have used the values Rin = 2i?max — 5 
and Rout = 2i?max, where i^max is the maximum radius 
of the cluster. The random walk experienced by the in- 
coming particles is affected by the interactions exerted 
by the particles already attached to the cluster. We have 
taken this fact into account by using a Metropolis algo- 
rithm inspired by Refs. Suppose that the cluster 
is composed by n — 1 particles, placed at the points r^, 
i = l,...,n— 1. At some time t the incoming particle 
occupies the position fn and has an interacting energy 
E = Yl^=i ^ni • At time t + 1 we compute a new position 
r^; the particle arrives there by means of jump of length 
(i in a direction chosen at random. The movement to is 
performed rigidly, without changing Un. The energy ex- 
perienced in the new position is , and the total change 
in the energy due to the movement is = E^ — E. If 
< 0, then the movement is accepted; if AE > 0, we 
compute the quantity p = exp(— A£^/Tr), where 



(2) 



{Tr is a reduced temperature^ related to the intensity of 
the interaction and the actual temperature T.) In this 
latter case, the movement is accepted with probability 
p. After every accepted movement, the moment of the 
random walker is oriented along the direction of the total 
field on its position. This fact indeed assumes that the re- 
laxation time for the orientation of particles is very short 
in comparison with the movement of the center of mass. 
The particle sticks to the cluster when it overlaps one or 
more particles already incorporated. After sticking, the 
newly attached particle undergoes one last relaxation. 



III. FRACTAL PROPERTIES OF THE CLUSTERS 

The purpose of this Section is to analyze the fractal 
properties of the clusters generated using the prescrip- 
tion outlined above. In Fig.|l| we have represented typical 
clusters of 1000 particles, grown at four different values 
of Tr. The effects of temperature can be seen by compar- 
ing with a pure non-interacting DLA cluster, as shown 
in FigJ^(d). At low T^, Fig.|l|(a), the dipolar clusters 
have a lesser branched and more open structure, that 
is to say, they have a smaller ratio of bifurcation. Even 
though the clusters may seem to be anisotropic, they still 
posses spherical symmetry: When collapsing an ensem- 
ble of clusters at the same Tp, we recover a perfectly 
symmetric structure. When increasing the temperature, 
the branching of the clusters increases correspondingly; 



for Tr = 10 (the largest value simulated), Fig.|i|(c), the 
clusters are completely indistinguishable from true DLA. 
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FIG. 1. Typical 2D dipolar clusters of 1000 particles, gen- 
erated for several values of Tr. (a) Tr — 0^ D — 1.20 =b 0.02. 

(b) Tr = 10"^ D = 1.35±0.04. (c) Tr = 10, D = 1.74±0.02. 

(c) Pure DLA cluster, grown in the limit Tr cxo, 
D = 1.71 ±0.01. 

In order to quantify this temperature dependence, we 
have computed the fractal dimension D of the clusters, 
determined from a log-log plot of the radius of gyration 
as a function of the number of particles N pi. 
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The algorithm was checked by computing the dimension 
of an ensemble of 100 clusters of 1000 particles, generated 
in the limit ^ oo (pure DLA). The value computed 
was D = 1.71 ± 0.01; we recover, within the error bars, 
the well-known result D = 1.715 ± 0.004 

Fig.^ shows a plot of I) as a function of Tr for the whole 
range of values analyzed. It seems to exhibit a smooth 
increasing behavior when increasing T^, the values of D 
ranging between the limits corresponding to = (ag- 
gregation with dipolar forces of infinite strength or zero 
temperature) with D = 1.13±0.01 and = oo (aggre- 
gation with no interactions or infinite temperature), with 
D = 1.71 ± 0.01. We can compare our results with the 
colloidal aggregation experiment described in |l^ , where 
magnetic particles were employed for which T^ ^ '^^^ 
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at room temperature. In our simulations, the fractal 
dimension obtained for clusters grown at Tr = 10~^ is 
D = 1.35 ± 0.04, a value clearly different from the one 
computed at = 0. 
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FIG. 2. Fractal dimension as a function of Tr for dipolar 
clusters. 

The shape of Fig.| can be explained as follows. Given 
the expression (|l|) for the dipolar energy, it is energet- 
ically favorable for an incoming particle to stick at the 
tip of a branch, contributing thus to its growth, rather 
than sticking on any side and splitting it. At low tem- 
peratures, therefore, the most probable scenario is that 
of a cluster with a very small branching ratio and a low 
fractal dimension. This in consistent with the fact that, 
at = 0, the fractal dimension seems to decrease when 
increasing N; we cannot even reject the possibility of 
a dimension equal to 1 in the limit ^ oo. On the 
other hand, at large temperature the growth of a given 
branch and its split are equally likely events. The effect 
of dipolar interactions is overcome by thermal disorder 
and we recover the original DLA model with no inter- 
actions. At intermediate values of T^, there would be a 
competition between growth and splitting, ruled by the 
thermal disorder as well as the dipolar interactions. We 
should then expect a continuous change in the geometry 
of the cluster and, therefore, a smooth dependence of D 
on the temperature. 

The fractal dimension we have considered so far is 
a macroscopic property of dipolar aggregates. How- 
ever, their composing particles also bear a rigid magnetic 
dipole, which confers a microscopic structure to the clus- 
ters. In order to obtain information about this structure, 
we have analyzed the relative orientation between pairs 
of dipoles. To this end, we define the function gri^) as 
the probability density that the relative angle formed by 
the directions of a pair of dipoles randomly chosen from 
a cluster at temperature Tr is included in the interval 

-\- dO]; is defined to be normalized to the interval 
[0,7r]. In practice, if n(T^, 6) is the number of pairs with 



a relative angle between and 0-\-A0, for a fixed angular 
increment A^, then we have 

N being the number of particles in the cluster. gT is 
a measure of the order of the dipoles on the cluster. 
In a completely disordered distribution, all relative ori- 
entations are equally probable and, therefore, we have 
g{6) = I/tt. On the other hand, in a distribution in 
which all the dipoles point in the same direction (for 
example, in the presence of a strong magnetic field) ev- 
ery pair forms a relative angle of zero radians, and thus 
g{0) = S{0). 

Fig.^ depicts gri^) computed from several ensembles 
of clusters grown at different values of Tr , between and 
10. 
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FIG. 3. Orientational correlation griO) as a function of 
cos^, computed for dipolar clusters at different values of Tr. 
At high temperatures the slope is almost zero. Solid lines are 
least-squares fittings. 

Numerically we observe that ^t(^) fits extremely well 
to the function 

gT{0) =a^b{Tr)cosO. (5) 

The normalization condition of gr implies obviously a = 
I/tt, as observed. On the other hand, b{Tr) seems to be 
a decreasing function of temperature, between the limits 
6(0) = 0.0594 ± 0.0005 and b{oo) - 0. The slope b{Tr) 
can be seen as a measure of the degree of internal order 
in the orientation of the dipoles in the cluster. In the low 
temperature limit (in the absence of any thermal disor- 
der) the orientation of the dipoles is strongly correlated, 
and this fact is reflected in a nonzero slope b. As tem- 
perature raises, griO) decreases continuously. In the high 
temperature limit dipoles show a large disorder, imposed 
by the intrinsic fractal geometry, and the function gr is 
almost flat (all relative orientations are equally proba- 
ble). Once again, the apparently continuous variation of 
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gT hints towards a smooth dependence of the geometry 
of the clusters on the temperature. 



IV. CONCLUSIONS 

We have investigated the effects of dipolar interactions 
in particle-cluster aggregation in two dimensions. The 
relevant parameter in the model is the dimensionless tem- 
perature Tr defined in which relates the real temper- 
ature T and the strength of the magnetic interactions 
/i. At low temperatures we observe clusters with a less 
branched and more open structure than free DLA, and 
correspondingly, a fractal dimension close to 1. This is 
an effect of the dipolar interaction on the local growth- 
site probability distribution [||. The growth of a given 
branch is a more likely event than its split, resulting in 
an enhancement of the screening of the inner regions. At 
high temperatures the fractal dimension raises its value 
until it reaches the limit of free DLA. The plot of I) as a 
function of T^, FigJ|, seems to show a smooth behaviour, 
resulting from the competition between dipolar attractive 
forces and thermal disorder. The internal structure of the 
clusters, probed through the function ^t, seems also to 
show a smooth transition between an ordered state at 
low temperature, with long range correlations between 
dipoles, and a disordered state at high temperature, in 
which all relative orientations are equally probable. This 
loss of order is explained as the effect of the geometrical 
disorder induced by the fractal character of the clusters. 
Our results could be relevant to understand the processes 
of cluster aggregation in dipolar colloids. Even though 
our result are relative to attractive dipolar interactions, 
they could also be applied to study cluster growth in 
Langmuir monolayers; the similarity with experimental 
clusters has already been pointed out in Ref. [Q, in the 
particular case of isotropic dipolar interactions (potential 
decaying with distance as r~^). 

As a final remark, we would point out that, even 
though all our results are consistent with a geometry 
continuously varying with the temperature, they are also 
compatible with a different scenario, in which there is 
sharp crossover between the low temperature dimen- 
sion 1^0 — 1-13 and the high temperature dimension 
Doo — 1.71. From Fig.^, the crossover temperature could 
be estimated to be lO"^-^ - 3 x 10"^. The pres- 

ence of the sharp fall in D would then be smoothed out 
by finite-size effects, unavoidable for the cluster sizes we 
are considering. Further work should be done in order to 
elucidate this possibility, especially by simulating aggre- 
gates larger than those actually available with our present 
computer resources. 
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